Ballistic Thermal Conduction across Acoustically Mismatched Solid Junctions 
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We derive expressions for energy flow in terms of lattice normal mode coordinates and energy 
transmission involving reduced group velocities. With a version of Landauer formula appropriate 
for lattice dynamic approach, the phonon transmission coefficients and thermal conductance are 
calculated for two kinds of acoustically mismatched junctions: different chirality nanotubes (11,0) 
to (8, 0), and Si-Ge superlattice structure. Our calculation shows a mode-dependent transmission in 
nanotube junction and a resonantly modulated ballistic thermal conductance in superlattice. The 
superlattice result suggests a new interpretation of the experimental data. Our approach provides 
an atomistic way for the calculation of thermal conduction in nanostructure. 
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Rapid progress in the synthesis and processing of ma- 
terials with structures of nanometer length scales has cre- 
ated a demand for the understanding of thermal trans- 
port in nano-scale low dimensional devices 0, 0, S • 
Nanostructures offer a new way of controlling thermal 
transport by tuning dispersion relations and other pa- 
rameters 0. Recent experimental and theoretical stud- 
ies have revealed novel features of phonon transport in 
these systems, such as the size-dependent anomalous heat 
conduction in one-dimensional (ID) chain 2] and the uni- 
versal quantum thermal conductance [j| • Thermal trans- 
port in nanostructures may differ from the predictions 
of Fourier's law based on bulk materials; this may hap- 
pen because of the existence of many acoustically mis- 
matched interfaces in nanostructures and because the 
phonon mean free path is comparable to the size of the 
structure 0. An understanding of the thermal conduc- 
tion across acoustically mismatched solid interfaces is a 
necessary requirement for thermal transport engineering. 

The study of thermal transport across interfaces dates 
back as early as to 1940s when Kapitza resistance (£( 
was reported and much work has been done in this field 
4, 5j. In general, theoretical modeling of this prob- 
lem has been undertaken either by the acoustic-mismatch 
model (AM) with scalar elastic waves, or by the diffuse- 
mismatch model (DM) with Boltzmann transport equa- 
tion 6]. Some numerical methods such as molecular dy- 
namic simulation pfl have also been used. While AM and 
DM models provide some useful reference calculations, 
scalar wave model and Boltzmann transport equation 
are only phenomenological descriptions and they have 
ignored the complexity of the interface. Atomic-level lat- 
tice dynamic (LD) approach should be the right way of 
capturing the mechanism of heat transport. However, 
after its early proposal in Ref. |8|, there has been little 
further work using this approach. 

Landauer formula for ballistic heat transport has been 
used || for the prediction of universal quantum heat con- 
ductance at very low temperatures. The formula derived 
under continuum assumption cannot be applied straight- 



forwardly to systems on nanoscale where atomic details 
are important. In this paper, we outline a new derivation 
from the point of view of lattice dynamics. The lattice 
formulation takes into account the different masses of 
the atoms and various vibrational modes. We also pro- 
pose a method of computing the transmission coefficients 
by solving a set of dynamical equations with scattering 
boundary conditions. The method is applied to compute 
the transmission coefficients of a carbon-nanotube junc- 
tion and superlattice structures, and the thermal conduc- 
tance is obtained with Landauer formula. 

We consider systems with perfect leads on the left and 
right with an arbitrary interaction at the junction. The 
Hamiltonian for such a system of vibrating atoms under 
linear approximation is given by 
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where I or V denotes a unit cell, i,j the position in a 
cell, a, /3 the direction of vibrating motions of atoms, and 
uf i the displacement from equilibrium of the atom 
with equilibrium position Rj,j. A local energy density 
can be defined through the energy in cell I, as p(r) = 
£iS(r— Rj) where R; is a lattice vector. An expression 
for the heat current in the z direction can be derived from 
the energy continuity equation, dp/dt + V- j = 0, as 
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In this equation, q z , f2(q, t) and L z are the z-component 
momentum, the energy density in momentum space and 
the system length along z direction, respectively, with 
time average represented by ( ) t . The energy density in 
momentum space is given by 
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where Qq's are normal mode coordinates for 
the n-th branch phonons, with mode uf in (t) — 

S= g ""( ( l) ei(q ' R, " Ut) ' Where e«„(q)e7«,(q) = 5„,n' 

[lOj. Combining Eq. (J2J and Eq. J3J, the time-averaged 
energy current along z-direction is 
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The transmission coefficients arc defined with respect to 
the normalizations of incoming and outgoing waves. The 
total energy current from one particular lead, say the 
left lead, is an arbitrary superposition of all the modes. 
Thus, the motion of the atoms is described by the wave: 



+ E( E <,^,W;f>,-q'), (5) 
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where is the amplitude of the mode n in lead <r, while 
tnn' ls transmission/reflection amplitude from mode 
n' in lead a' to mode n in lead a. Note that q and q' 
satisfies u> = w n (q) = w n (q'). Similar expression can be 
written down for the right lead. Since the energy is con- 
served, there is no net energy accumulation in the junc- 
tion, which means that the time-averaged energy currents 
from both sides are equal, I L = I R . This condition leads 
to the following identity for the transmission amplitudes: 
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The important difference with the continuum case |9| is 
that we need to replace the group velocity with a re- 
duced group velocity v° = v^/F z , where F z is the length 
of unit cell along the z direction, — du°/dq z . For ID 
quantum thermal energy current, with the help of the 
relation (JBJ, after quantizing Eq. we get the formula 
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v„/Vm is the energy transmis- 



sion probability, f(u>,T a ) is the Bose- Einstein distribu- 
tion for the left or right lead. Eq. Q is the Landauer 
formula for quantum thermal energy flow. The thermal 
conductance is obtained by taking an infinitesimal tem- 
perature difference. 

The central issue now is to have an efficient method to 
calculate the transmission coefficients across junctions at 
the atomistic level. This appears to be a difficult problem 
8] in general taking into account the complexity of the 
interface. Transfer matrix method can be used for simple 
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FIG. 1: (a) Structure of a (11,0) and (8,0) nanotube junc- 
tion, (b) The energy transition T n (uj) as a function of angular 
frequency u. The inset shows the dispersion relation of the 
infinite-long (11,0) nanotube. 



ID models with a few interfaces. But this method could 
not be applied to a large ID system or 3D systems due 
to numerical instability. 

We propose what we call the scattering boundary 
equation (SBE) method as a numerically more satisfac- 
tory solution. Each atom in the system satisfies the 
dynamic equation, — rrii^Ui + J^j Kij( u i ~ u j) = 0, 
where the 3x3 matrix Kij is the force constants be- 
tween atom i and j. The boundary conditions are of 
the form u 9i n' + ^2 n r n u_ ? ,„ for the incoming waves and 
J2 n *n"g',n f° r tnc outgoing waves, where u 9i „ and \i' q , n 
are the eigen modes on the left and right leads, while 
the reflection and transmission coefficients r n and t n are 
unknown. These equations in matrix form are illustrated 
as 
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where Ka = —niiUj + ^2jKij. In ID, the matrix is 
square and can be both analytically and numerically 
solved by the conventional method. For higher dimen- 
sions, however, the boundary conditions are complicated 
and the number of equations may be larger than that of 
variables. But these equations are not linearly indepen- 
dent and can still be solved by the singular value decom- 
position method. 

We first discuss the SBE results of the transmission co- 
efficients for a nanotube junction, shown in Fig. ^ This 
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calculation involves 504 equations and 455 variables. The 
semiconductor nanotube junction structure is optimized 
by a second-generation Brenner potential with the 
force constants derived from the same potential. The 
phonon dispersion for nanotube (11,0) is illustrated in 
the inset, in which four acoustic branches are considered: 
the longitudinal mode (LA), doubly degenerate trans- 
verse mode (TA), and the unique twist mode (TW) in 
nanotubes. Although all modes of a given frequency are 
considered, we did not find mode-mixing behavior among 
acoustic modes at the lower frequency range. The trans- 
mission for LA mode stays around 0.8 with only small 
changes. This value is below the AM model prediction 
of 0.98 with the longitudinal group velocity 20.18 km-s -1 
and 20.95 km-s" 1 for (11,0) and (8,0). In contrast, the 
transmission for TA mode decreases with frequency. This 
can be accounted for by a nearly quadratic dispersion re- 
lation of the TA mode, as illustrated in the inset. The 
transmissions of the TW mode and many other optical 
modes are nearly zero or very small. This appears related 
to the difference in rotational symmetries of these modes. 
We propose that this kind of mode-dependent transmis- 
sion behavior may be important for further application 
such as phonon filters. 

Next, we consider two acoustically mismatched chains 
of the same atomic mass connected with springs of differ- 
ent stiffnesses K\ , and lattice constants a, b on each 
side. The energy transmission by transfer matrix calcu- 

lation is \t\ 2 U = 4 (§* ) 2 /(l + ff)' where v L = v L /a, 

vr = vn/b. Using mass densities p^ — m/a, pa = m/b 
and zl — plvl, zr — prvr, the AM model gives |l2( 

-, iz ^ ZR y> = 4 (SA-) l(l + i n agreement with our 

(zl+zrY \ Vl J \ VL J 

result. Thus, for simple lattice structure and smooth 
junction, the results of LD and AM are the same. 

For a more realistic model, not only the physical lat- 
tice may be composite, but also the boundaries can be 
more complex. We choose effective stiffness to fit the dia- 
mond structure Silicon (Si) and Germanium (Ge) phonon 
dispersion relations. For Si and Ge along the direction 
T-X, the highest frequency for acoustic phonon branch 
is about 400 cm -1 and 250 cm" 1 , respectively |l3j]. If 
the phonon dispersion relation along this direction can 
be regarded as one-dimensional 0], the spring stiffness 
for each can be calculated by w — ^Ak/m = \f\K. We 
use these experimental parameters K2 = 2.QK\ for Si 
and Ge in the following discussion. 

Fig-Elshows how the energy transmission changes with 
the boundary spring stiffness Kf. It is reasonable to as- 
sume that the stiffness Kf for the Si-Ge junction is be- 
tween pure Si and Ge lattices. When 7 = Kf/K\ = 1.5 
or 2.3, the transmission is higher than the prediction 
given by direct acoustic mismatch model. Moreover, 
when the lattice structure is not a simple lattice, we also 
find that AM model does not hold any more. The in- 
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FIG. 2: Energy transmission T(ui) as a function of scaled 
angular frequency lo/ \f 4 A'i . Different curves correspond to 
different coupling strength ratios 7 = Kf/Ki at the junction 
with the spring stiffness parameters K2 = 2.6-Ki. The inset 
shows the difference in the prediction of energy transmission 
for the composite lattice structure, K2 = 0.3846 A'i, K3 — K\ 
for AM-1 and LD-1; K 2 = 1.1538-ffi, K 3 = 3.0iG for AM-2 
and LD-2. 

set in Fig. [21 shows the discrepancy in the prediction of 
the energy transmission by AM and LD model for a sim- 
ple lattice with stiffness K\ on the left and alternating 
stiffness K2 and K3 on the right. Comparing with LD 
result, the prediction of AM model overestimates the en- 
ergy transmission at high frequency for mismatched com- 
posite lattice structure. 

Furthermore, when the interfaces are made of period- 
ical superlattices of Si and Ge, interesting phenomenon 
emerges. Fig. ^a) shows a series of peaks appearing 
in the transmission for the superlattice of alternating Si 
and Ge monolayer. We find that this behavior can be 
understood as a resonance effect. This kind of resonant 
tunneling transmission has been reported in |l4j . When 
the superlattice layer is composed of several monolayers 
of Si and Ge, normal mode calculation shows that reso- 
nant frequencies fall into different parts and band gaps 
are brought about by the repetition of superlattice struc- 
ture. So it can be seen that the total energy transmission 
is the acoustic mismatch transmission modulated by the 
resonance. 

We find that these resonance gaps will relate the ther- 
mal conductance of superlattice with the thickness of the 
period layer. We map the ID chain results to a 3D lat- 
tice using an average lattice constant Iq — 5.5 A for the 
cross section of the cell, and 8 chains per conventional 
diamond-type cell. Our result of thermal conductance 
is given as G s» 10 9 Wm" 2 K" 1 at 200 K, which a gree s 
quantitatively with the measured effective values in [15j . 
From Fig. we see that thermal conductance decreases 
with the increase in the thickness of period L at high tem- 
peratures (> 100 K) due to the appearance of mini gaps 
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FIG. 3: Energy transmission T(u) of a superlattice with the 
Si and Ge leads; u> is scaled by y/4Ki. (a) Superlattice of 
alternating Si and Ge monolayer, n is the number of periods, 
(b) Superlattice of ten period layers, each layer with 5 x 5 Si 
and Ge monolayers. 
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FIG. 4: Thermal conductance as a function of temperature. 
Curves labeled 10 and 40 represent ten superlattice period 
layers, each layer with 5x5 and 20 x 20 Si and Ge monolayers, 
respectively; For comparison, the acoustic mismatch (AM) 
model result of a single Si-Ge interface is also plotted. 



in the transmission, which agrees with the experimen- 
tal result in 0. We argue that in the several superlat- 
tice length scale the thermal conductance is wave-natured 
ballistic transport because the phonon mean free path is 
comparable to this len gth scale. But for the samples of 
thickness of order /im |l5j, the total thermal resistance 
will be the diffusive summation of the microscopic bal- 
listic thermal resistance because wave coherence will be 



destroyed after transporting across a few superlattice pe- 
riods. So the effective thermal conductance in will 
exhibit microscopic ballistic thermal transport behavior. 
Authors of Ref [16| interpret this kind of thermal con- 
ductance as "partially coherent heat conduction" by in- 
troducing an imaginary wave vector. When the thickness 
of the superlattice is larger than the phonon mean free 
path, the conventional diffusive conductance comes into 
play and the thermal conductivity will rise with increas- 
ing L 15]. These qualitative agreements between our 
calculations and the experimental results suggest a bal- 
listic heat transport as the dominating mechanism within 
a few superlattice periods. Another feature we find from 
our calculation is that for lower temperatures (< 50 K), 
the superlattice conductance becomes independent of the 
thickness and tends to some universal value. We expect 
that some experiments will verify this point in future. 

In summary, we have derived a Landauer formula suit- 
able for lattice dynamic calculations. The transmission 
coefficients are calculated by a linear equation solver with 
proper boundary conditions. We have applied our meth- 
ods to two types of junction structures, and give signifi- 
cant results. 
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